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Abstract 

The prompt production of the charmonium Xci and Xc2 mesons has been studied 
in proton-proton collisions at the Large Hadron Collider at a centre-of-mass energy 
of -^/s = 7 TeV. The Xc mesons are identified through their decays Xc J/i'l 
with J/ if) — /i^ /U~ using 36 pb~^ of data collected by the LHCb detector in 
2010. The ratio of the prompt production cross-sections for the two Xc spin states, 
o'(Xc2) / o"(Xci), has been determined as a function of the J/tp transverse momentum, 
p'^^ , in the range from 2 to 15 GeV/c. The results are in agreement with the next- 
to-leading order non-relativistic QCD model at high p'!^^ and lie consistently above 
the pure leading-order colour singlet prediction. 
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1 Introduction 



Explaining heavy quarkonium production remains a challenging problem for Quantum 
Chromodynamics (QCD). At the energies of the proton-proton [pp) collisions at the Large 
Hadron Collider, cc pairs are expected to be produced predominantly via Leading Order 
(LO) gluon-gluon interactions, followed by the formation of the bound charmonium states. 
While the former can be calculated using perturbative QCD, the latter is described by 
non-perturbative models. Other, more recent, approaches make use of non-relativistic 
QCD factorization (NRQCD) which assumes a combination of the colour- singlet (CS) 
and colour-octet (CO) cc and soft gluon exchange for the production of the final bound 
state [T]. To describe previous experimental data, it was found to be necessary to in- 
clude Next-to-Leading Order (NLO) QCD corrections for the description of charmonium 
production [21 13]. 

The study of the production of P-wave charmonia Xcj{^P), with J = 0, 1, 2, is impor- 
tant, since these resonances give substantial feed-down contributions to the prompt J/ip 
production through their radiative decays Xc ~^ J/i^ 7 and can have significant impact 
on the measurement of the J/ip polarisation. Furthermore, the ratio of the production 
rate of Xc2 to that of Xci is interesting because it is sensitive to the CS and CO production 
mechanisms. 

Measurements of Xc production and the relative amounts of the Xci and Xc2 spin states, 
have previously been made using different particle beams and energies [U |5l [6] . In this 
Letter, we report a measurement from the LHCb experiment of the ratio of the prompt 
cross-sections for the two Xc spin states, o"(Xc2) / o'ixa), as a function of the J/ip transverse 
momentum in the range 2<p^^ < 15 GeV/c and in the rapidity range 2.0 <y-'^^ <4.5. 
The Xc candidates are reconstructed through their radiative decay Xc — ^ J/'ip 7, with 
J/V' — )■ using a data sample with an integrated luminosity of 36 pb~^ collected 

during 2010. In this Letter, prompt production of Xc refers to Xc mesons that are produced 
at the interaction point and do not arise from the decay of a 6-hadron. The sample 
therefore includes Xc from the decay of short-lived resonances, such as ip{2S), which are 
also produced at the interaction point. All three Xcj states are considered in the analysis. 
Since the Xco J/'^P 7 branching fraction is ~ 30 times smaller than those of the Xci and 
Xc2, the yield of Xco is not significant. The measurements extend the p'!^' coverage with 
respect to previous experiments. 

2 LHCb detector and selection requirements 

The LHCb detector [7] is a single-arm forward spectrometer with an angular coverage from 
approximately 10 mrad to 300 mrad (250 mrad) in the bending (non-bending) plane. The 
detector consists of a vertex detector (VELO), a dipole magnet, a tracking system, two 
ring-imaging Cherenkov (RICH) detectors, a calorimeter system and a muon system. 

Of particular importance in this measurement are the calorimeter and muon systems. 
The calorimeter consists of a scintillating pad detector (SPD) and a pre-shower, followed 
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by electromagnetic (ECAL) and hadronic calorimeters. The SPD and pre-shower are 
designed to distinguish between signals from photons and electrons. The ECAL is con- 
structed from scintillating tiles interleaved with lead tiles. Muons are identified using hits 
in detectors interleaved with iron filters. 

The signal simulation sample used for this analysis was generated using the Pythia 6.4 
generator [8j configured with the parameters detailed in Ref. [9]. The EvtGen (10] . 
Photos pjj and Geant4 [12] packages were used to decay unstable particles, generate 
QED radiative corrections and simulate interactions in the detector, respectively. 

The trigger consists of a hardware stage followed by a software stage which applies a full 
event reconstruction. For this analysis the trigger selects a pair of oppositely charged muon 
candidates, where either one of the muons has a transverse momentum pt>1-8 GeV/c 
or one of the pair has > 0.56 GeV/c and the other has > 0.48 GeV/c. The invariant 
mass of the candidates is required to be greater than 2.9 GeV/c^. The photons are not 
involved in the trigger decision for this analysis. 

Photons are identified and reconstructed using the calorimeter and tracking systems. 
The identification algorithm provides an estimator for the hypothesis that a calorimeter 
cluster originates from a photon. This is a likelihood-based estimator constructed from 
variables that rely on calorimeter and tracking information. For example, in order to 
reduce the electron background, candidate photon clusters are required not to be matched 
to a track extrapolated into the calorimeter. For each photon candidate a likelihood (CL^) 
is calculated based on simulated signal and background samples. The photons identified 
by the calorimeter and used in this analysis can be classified as two types: those that have 
converted in the material after the dipole magnet and those that have not. Converted 
photons are identified as clusters in the ECAL with correlated activity in the SPD. In order 
to account for the different energy resolutions of the two types of photons, the analysis is 
performed separately for converted and non-converted photons and the results combined 
as described in Sect. [3j Photons that convert before the magnet require a different analysis 
strategy and are not considered here. The photons used to reconstruct the Xc candidates 
are required to have a transverse momentum > 650 MeV/ c, a momentum p'~' >5 GeV/ c 
and a likelihood CL^ > 0.5. 

The muon and J/ip identification criteria are identical to those used in Ref. ||13] : 
each track must be identified as a muon with > 700 MeV/ c and a quality of the track fit 
X^/ndf < 4, where ndf is the number of degrees of freedom. The two muons must originate 
from a common vertex with a probability of the vertex fit > 0.5%. In addition, in this 
analysis the invariant mass is required to be in the range 3062 — 3120 MeV/c^. The 

J/ip pseudo-decay time, tz, is used to reduce the contribution from non-prompt decays, 
by requiring = {zj/^ — zpy)Mj/^ / Pz< 0.1 ps, where Mj/^ is the reconstructed dimuon 
invariant mass, zj/^ — zpy is the z separation of the reconstructed production and decay 
vertices of the dimuon, and Pz is the ^-component of the dimuon momentum with the 
z-axis parallel to the beam line. Simulation studies show that, with this requirement 
applied, the remaining fraction of Xc from 6-hadron decays is about 0.1%. This introduces 
an uncertainty much smaller than any of the other systematic or statistical uncertainties 
evaluated in this analysis and is not considered further. 
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M(n* n" T) - M{n* [MeV/tf] M{n* \i y) - M(n* n") [MeV/c^] 



Figure 1: Distribution of AM = M(/i+/i j) — M {fi'^ ) for selected candidates with 
3<Pt^ < 15 GeV/c for (a) converted and (b) non-converted photons. The lower solid curves 
correspond to the XcO, Xci and Xc2 peaks from left to right, respectively (the XcO peak is barely 
visible). The background distribution is shown as a dashed curve. The upper solid curve corre- 
sponds to the overall fit function. 



In the data, the average Xc candidate multiphcity per selected event is 1.3 and the 
percentage of events with more than one genuine Xc candidate (composed of a unique 
J/ip and photon) is estimated to be 0.23% from the simulation. All Xc candidates are 
considered for further analysis. The mass difference, AM = M {^~^ /i~ 7) — M {fi~^ /i~), of 
the selected candidates is shown in Fig. [l] for the converted and non-converted samples; 
the overlaid fits are described in Sect. HI 



3 Experimental method 

The production cross-section ratio of the Xc2 and xa states is measured as 

^iM = ^ ^ J/i^ 7) . . 

where i3(Xci — ^ J/i^ l) and B{xc2 J/i^ l) are the Xci and Xc2 branching fractions to 
the final state J/ip 7, and 



where e^J^ is the efficiency to trigger, reconstruct and select a J/ip from a Xcj decay, e^"^ is 
the efficiency to reconstruct and select a photon from a Xcj decay and e^^^ is the efficiency 
to subsequently select the Xcj candidate. 

Since the mass difference between the xa and Xc2 states is 45.54 ±0.11 MeV/c^, the 
signal peaks cannot be separately isolated using the calorimeter information. An unbinned 
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maximum likelihood fit to the AM mass difference distribution is performed to obtain 
the three N^^j yields simultaneously. The determination of the efficiency terms in Eq. [2] 
is described in Sect. 13. 1[ 



The signal mass distribution is parametrised using three Gaussian functions (J^,'ig for 
J = 0, 1, 2). The combinatorial background is described by 

•^bgd = a;"(l-e'?(i-^')) (3) 

where x = AM / mo and mo, a, b and c are free parameters. 

A possible source of background from partially reconstructed decays is due to 
ip{2S) — > J/tp 7r° tt" decays where the J/ip and a photon from one of the neutral pions 
are reconstructed and selected as a Xc candidate. Simulation studies show that the ex- 
pected yield is ~ 0.1% of the signal yield and this background is therefore neglected for 
this analysis. 

The overall fit function is 



Slg 

J=0 



J=0 



^bgd, (4) 



where f^^j are the signal fractions. The mass differences between the Xci and Xc2 states 
and the Xci and Xco states are fixed to the values from Ref. [T^ . The mass resolutions for 
the Xc states, cries' , are given by the widths of the Gaussian functions for each state. The 
ratios of the mass resolutions, a^^^ / a^^^ and a^g° / a^^^ , are taken from simulation. The 
value of (T^f / cr^'gi is consistent with the value measured from data, fitting in a reduced 
AM range and with a simplified background parametrisation. 

With the mass differences and the ratio of the mass resolutions fixed, a fit is performed 
to the data in the range 3 <p^^ < 15 GeV/c, in order to determine the Xci mass resolution 
(J^es- This range is chosen because the background has a different shape in the p'^^ 
bin 2 — 3 GeV/ c and is not well described by J-bgd when combined with the rest of the 
sample. Simulation studies show that the signal parameters for the Xcj states in the p'^^ 
bin 2 — 3 GeV/c are consistent with the parameters in the rest of the sample. The AM 
distributions of the fits to the converted and non-converted candidates are shown in Fig. [T] 
The mass resolution, cx^g^, is measured to be 21.8 ± 0.8 MeV/c^ and 18.3 ± 0.4 MeV/c^ for 
converted and non-converted candidates respectively. The corresponding values in the 
simulation are 19.0 ±0.2 MeV/c^ and 17.5 ±0.1 MeV/c^ and show a weak dependence of 
(T^'gi on p'^^ which is accounted for in the systematic uncertainties. 

In order to measure the Xc yields, the fit is then performed in bins of p^^ in the range 
2 < p'^^ < 15 GeV/ c. For each p^^ bin, the mass differences, the ratio of the mass resolu- 
tions and 0"^^^ are fixed as described above. In total, there are eight free parameters for 
each fit in each bin in p^^ and the results are summarized in Table 1 The total observed 
yields of Xco, Xci and Xc2 are 820 ±650, 38 630 ±550 and 26 110 ±620, respectively, cal- 
culated from the signal fractions f^^j and the number of candidates in the sample. The 
raw Xc yields for converted and non-converted candidates are combined, corrected for 
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Table 1: Signal Xc yields from the fit in each p'^^ bin. 



p-!!^ (GeV/c) 


Converted photons 


Non-converted photons 


Xci yield 


Xc2 yield 


X^ 1 ndf 


Xci yield 


Xc2 yield 


/ ndf 


2-3 


3120 ±248 


2482 ± 301 


0.91 


4080 ± 246 


3927 ±280 


1.02 


3-4 


3462 ± 224 


3082 ± 249 


0.81 


4919 ±183 


3443 ± 207 


1.02 


4-5 


3235 ± 146 


1769 ±174 


1.03 


4497 ±134 


2718 ±143 


1.08 


5-6 


2476 ±110 


1443 ±121 


0.84 


3203 ±105 


1999 ±107 


1.45 


6-7 


1497 ±80 


736 ± 89 


1.05 


1946 ± 78 


1338 ±83 


0.78 


7-8 


933 ± 77 


658 ±86 


0.77 


1342 ± 59 


747 ± 60 


1.15 


8-9 


660 ± 47 


302 ±51 


0.90 


817 ±43 


395 ± 42 


0.78 


9-10 


451 ±34 


142 ±35 


0.82 


501 ±32 


256 ±31 


1.09 


10-11 


255 ±25 


86 ±26 


1.13 


317±26 


188 ±25 


0.85 


11-12 


129 ±28 


99 ±30 


0.87 


222 ±19 


103 ±18 


0.93 


12-13 


129 ±16 


46 ±15 


1.09 


154 ±15 


50 ±13 


0.98 


13-15 


127 ±18 


42 ±20 


0.91 


158 ±18 


63 ±17 


1.05 



efficiency (as described in Sect. 3.1 ) and the cross-section ratio is determined using Eq. [T] 



3.1 Efficiencies 

The efficiency ratios to reconstruct and select Xc candidates are obtained from simulation. 
Since the photon interaction with material is not part of the event generation procedure, 
the individual efficiencies for converted and non-converted candidates are not separated. 
Therefore, the combined efficiencies are calculated. The ratios of the overall efficiency 
for the detection of J/V' mesons originating from the decay of a Xci compared to a Xci-, 
^jfl) I 1 consistent with unity for all "pl^^ bins, as shown in Fig. |2| The ratios of 
the efficiencies for reconstructing and selecting photons from Xc decays and then selecting 
the Xci ^^''^e^ef / ^I'^^^se?' ^Iso shown in Fig. [2j In general these efficiency ratios are 
consistent with unity, except in the bins 2 — 3 GeV/ c and 3 — 4 GeV/ c where the 
reconstruction and detection efficiencies for Xc\ are smaller than for Xc2- The increase in 
the efficiency ratio in these bins arises because the photon px spectra are different for Xc\ 
and Xc2- The photon > 650 MeV/c requirement cuts harder in the case of the Xc\ and 
therefore lowers this efficiency. The increase in the efficiency ratio is a kinematic effect, 
rather than a reconstruction effect, and is well modelled by the simulation. 
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Figure 2: Reconstruction and selection efficiency ratios in bins of . The ratio of the J/ip effi- 
ciency (e j/^ / e j/^ ) is shown with red circles. The ratio of the photon reconstruction and selection 
efficiency times the Xc selection efficiency {^^'''^(^^\ / ^'^"^^^eX^ shown with blue triangles. 

3.2 Polarisation 

The production of polarised Xc states would modify the efficiencies calculated from the 
simulation, which assumes unpolarised Xc- A measurement of the Xc polarisation would 
require an angular analysis, which is not feasible with the present amount of data. Various 
polarisation scenarios are considered in Table [2j Assuming no azimuthal dependence in 
the production process, the Xc — ^ J/iJ 7 system is described by three angles: 6j/^ , 6-^^ 
and 0, where 9ji^ is the angle between the directions of the muon in the J/ip rest frame 
and the J/ip in the Xc rest frame, 9-^^ is the angle between the directions of the J/ip in the 
Xc rest frame and the Xc in the laboratory frame, and is the angle between the plane 
formed from the Xc and J/ip momentum vectors in the laboratory frame and the J/ifj 
decay plane in the J/ip rest frame. For each simulated event in the unpolarised sample, 
a weight is calculated from the distribution of these angles in the various polarisation 
hypotheses compared to the unpolarised distribution. The weights in Table [2] are then 
the average of these per-event weights in the simulated sample. For a given (m^^j, m^^j) 
polarisation combination, the central value of the determined cross-section ratio in each 
p'^P' bin should be multiplied by the number in the table. The maximum effect from the 
possible polarisation of the xa and Xc2 mesons is given separately from the systematic 
uncertainties in Table |4] and Fig. [3| 
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Table 2: Polarisation weights inp^^ bins for different combinations of Xci and Xc2 polarisations, 
where rn^^j is the positive z-component of the total angular momentum J and Unpol. means 
the Xc is unpolarised. 



V A.CI 1 A.CZ / 


p^/^ (GeV/c) 


2-3 


3-4 


4-5 


5-6 


6-7 


7-8 


8-9 


9-10 


10-11 


11-12 


12-13 


13-15 


(Unpol,0) 


0.99 


0.97 


0.94 


0.91 


0.88 


0.87 


0.86 


0.86 


0.86 


0.85 


0.85 


0.88 


(Unpol,!) 


0.97 


0.98 


0.97 


0.95 


0.94 


0.94 


0.93 


0.93 


0.93 


0.93 


0.93 


0.93 


(Unpol,2) 


1.03 


1.04 


1.07 


1.11 


1.14 


1.17 


1.18 


1.18 


1.19 


1.18 


1.19 


1.16 


(0,Unpol) 


1.01 


0.99 


0.97 


0.93 


0.90 


0.89 


0.87 


0.86 


0.85 


0.87 


0.86 


0.84 


(l,Unpol) 


0.99 


1.00 


1.02 


1.04 


1.05 


1.06 


1.06 


1.07 


1.08 


1.07 


1.07 


1.08 


(0,0) 


1.00 


0.97 


0.91 


0.84 


0.80 


0.77 


0.75 


0.74 


0.72 


0.74 


0.74 


0.74 


(0,1) 


0.98 


0.97 


0.93 


0.88 


0.85 


0.83 


0.81 


0.81 


0.79 


0.81 


0.81 


0.78 


(0,2) 


1.04 


1.03 


1.03 


1.03 


1.03 


1.03 


1.03 


1.02 


1.00 


1.03 


1.03 


0.98 


(1,0) 


0.99 


0.97 


0.96 


0.94 


0.93 


0.92 


0.92 


0.92 


0.92 


0.91 


0.91 


0.95 


(1,1) 


0.97 


0.98 


0.98 


0.99 


0.99 


0.99 


0.99 


1.00 


1.00 


1.00 


1.00 


1.01 


(1,2) 


1.03 


1.04 


1.09 


1.15 


1.20 


1.23 


1.26 


1.26 


1.28 


1.26 


1.27 


1.25 



4 Systematic uncertainties 

The systematic uncertainties detailed below are measured by repeatedly sampling from 
the distribution of the parameter under consideration. For each sampled value, the cross- 
section ratio is calculated and the 68.3% probability interval is determined from the 
resulting distribution. The positive (negative) systematic uncertainty is defined as the 
difference between the upper (lower) bound of the interval and the central value of the 
cross-section ratio calculated using the procedure described in Sect. [3j 

The branching fractions used in the analysis are B{xci J/'ip 7) = 0.344 ± 0.015 and 
B{Xc2 J/'?/' 7) = 0.195 ± 0.008, taken from Ref. [13]. The systematic uncertainty on 
the cross-section ratio resulting from the Xc J/4' 7 branching fractions is determined 
by sampling from two Gaussian distributions, one for each of the xa and Xc2 branching 
fractions. The relative uncer_tainty is found to be in the range (5.2 — 7.4)% and is given 
for each bin of p^^ in Table 
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The simulation sample used to calculate the efficiencies has approximately the same 
number of Xc candidates as are observed in the data. The statistical errors from the finite 
number of simulated events are included as a systematic uncertainty in the final results. 
The uncertainty associated to this is determined by sampling the efficiencies used in Eq. [l] 
according to their errors. The relative systematic uncertainty due to the limited size of 
the simulation sample is found to be in the range (0.6 — 7.2)% and is given for each p'Jp' 
bin in Table [H 

The measured Xc yields depend on the values of the fixed parameters and the ffi range 
used. The associated systematic uncertainty has been evaluated by repeating the fit many 
times, changing the values of the fixed parameters and the fit range. Since the uncertain- 
ties arising from the fixed parameters are expected to be correlated, a single procedure 
is used simultaneously varying all these parameters. The systematic uncertainty is then 
extracted from the distribution of the Xc2 / Xci cross-section ratios given by these ffis as 
for the other systematic uncertainties. The Xc mass difference parameters are sampled 
from two Gaussian distributions with widths taken from the errors on the masses given 
in Ref. [Hj. The mass resolution ratios, cr^gf / a}^^ and a}^^ / cr^sS are varied according 
to the error matrix of the fit to the simulated sample in the range 3 <p'^^ < 15 GeV/c. 

The mass resolution a^^^ is also determined using a simplified background model and 
fitting in a reduced range. Simulation studies show that the value of a^^^ also has a weak 
dependence on p'Jp' . The value of a^^^ is randomly chosen to be either sampled from the 
default fit (described in Sect. [3]), according to its error, or sampled from the simplified fit, 
again according to its error, or by modifying it in each p'!p' bin according to the variation 
observed in the simulation. 

The systematic uncertainty associated with the shape of the ffited background function 
is incorporated by including or excluding the Xcq signal shape, which peaks in the region 
where the background shape is most sensitive. 

The background shape is also sensitive to the rise in the AM distribution. The 
systematic uncertainty from this is included by varying the lower edge of the fit range in 
the interval ± 10 MeV/ around its nominal value for each bin in p'^rp' . 

The overall systematic uncertainty from the ffi is then determined from the distribution 
of the cross-section ratio from repeating the sampling procedure described above many 
times. The relative uncertainty is found to be in the range (2.2 — 14.6)% and is given for 
each bin of p^^ in Table [s 

A systematic uncertainty related to the calibration of the simulation is evaluated 
by performing the analysis on simulated events and comparing the efficiency-corrected 
ratio of yields, (A^xc2 / ^xa) ' i^^"^ I ^'^''^)) the true ratio generated in the sample. A 
deviation of —9.6% is observed, caused by non-Gaussian signal shapes in the simulation 
from the calorimeter calibration. These are not seen in the data, which is well described 
by Gaussian signal shapes. No evidence for a deviation was observed when generating 
and ffiting distributions similar to those observed in the data, i.e. three overlapping 
Gaussians and a background shape similar to that in Fig. [Tj This deviation is included 
as a systematic error, by sampling from the negative half of a Gaussian with zero mean 
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Table 3: Summary of the systematic uncertainties on (t(xc2) / o'(Xci) in each pr^ bin. 



p^/^(GeV/c) 


2-3 


3-4 


4-5 


5-6 


6-7 


7-8 


Branching fractions 
Size of simulation sample 
Fit model 
Simulation calibration 


+0.08 
-0.09 

+0.01 
-0.01 

+0.04 
-0.05 

+0.00 
-0.08 


+0.08 
-0.08 

+0.01 
-0.01 

+0.05 
-0.04 

+0.00 
-0.07 


+0.06 
-0.06 

+0.01 
-0.01 

+0.03 
-0.03 

+0.00 
-0.05 


+0.07 
-0.07 

+0.01 
-0.01 

+0.03 
-0.03 

+0.00 
-0.05 


+0.07 
-0.06 

+0.02 
-0.01 

+0.03 
-0.04 

+0.00 
-0.06 


+0.07 
-0.06 

+0.02 
-0.02 

+0.05 
-0.04 

+0.00 
-0.06 




8-9 


9-10 


10-11 


11-12 


12-13 


13-15 


Branching fractions 
Size of simulation sample 
Fit model 
Simulation calibration 


+0.06 
-0.05 

+0.02 
-0.02 

+0.03 
-0.04 

+0.00 
-0.04 


+0.05 
-0.04 

+0.02 
-0.02 

+0.03 
-0.03 

+0.00 
-0.04 


+0.06 
-0.05 

+0.04 
-0.04 

+0.03 
-0.03 

+0.00 
-0.05 


+0.07 
-0.06 

+0.06 
-0.06 

+0.02 
-0.13 

+0.00 
-0.06 


+0.04 
-0.04 

+0.05 
-0.05 

+0.02 
-0.02 

+0.00 
-0.04 


+0.04 
-0.04 

+0.05 
-0.05 

+0.08 
-0.03 

+0.00 
-0.03 



and a width of —9.6%. The relative uncertainty on the cross-section ratio is found to be 
less than 6.0% and is given for each bin of p'!p' in Table 3 Other systematic uncertainties 
due to the modelling of the detector in the simulation are negligible. 

In summary, the overall systematic uncertainty is evaluated by simultaneously sam- 
pling the deviation of the cross-section ratio from the central value, using the distributions 
of the cross-section ratios described above. The systematic uncertainty is then determined 
from the resulting distribution as described in the introduction to Sect. |4j The separate 
systematic uncertainties are shown in bins of p'^^ in Table [3 and the combined uncertain- 
ties are shown in Table HI 



5 Results and conclusions 

The cross-section ratio, cr{xc2) / c"(Xci), measured in bins of p'!^^ is given in Table [i] and 
shown in Fig. [3] Previous measurements from WAll in 7r~Be collisions at 185 GeV/c gave 
o-(Xc2) / o-(Xci) = 1-4 ±0.6 S3], and from HERA-B in pk collisions at a/s = 41.6 GeV with 
p'J^'^ below roughly 5 GeV/c gave ct{xc2) I cr(Xci) = 1-75 ± 0.7 [5J. The data points from 
CDF [6j at ^/s = 1.96 TeV in pp collisions are also shown in Fig. |3} 

Comparisons to the theory predictions, in the rapidity range 2.0 < y'^^'^ < 4.5, from the 
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Table 4: Ratio a{xc2) / <^{Xci) in bins of prp' in the range 2<pr^ < 15 GeV/c and in the 
rapidity range 2.0 < 7/-^/^ <4.5. The first error is the statistical error, the second is the systematic 
uncertainty (apart from the branching faction and polarisation) and the third is due to the Xc — ^ 
J/ip 7 branching fractions. Also given is the maximum effect of the unknown Xc polarisations 
on the result as described in Sect. 13.21 



p-!!^ (GeV/c) 




Polarisation variation 


2-3 


1 oQ+0.12 +0.06 +0.08 
-'-•'J-^-0.13 -0.09 -0.09 


+0.06 
-0.05 


3-4 


1 oo+O.lO +0.03 +0.08 
-'-•'J^-0.09 -0.09 -0.08 


+0.06 
-0.05 


4-5 


1 CtO+^-Qf +0.04 +0.06 
^■^^-Om -0.06 -0.06 


+0.09 
-0.09 


5-6 


1 nS+'^-O'^ +0.04 +0.07 
-'-•Uo_o.o6 -0.06 -0.07 


+0.16 
-0.17 


6-7 


1 nO+'^-OS +0.03 +0.07 
-'-•UJ_o.09 -0.07 -0.06 


+0.22 
-0.22 


7-8 


1 nS+O-l^ +0.05 +0.07 
-'-•Uo_o.lO -0.07 -0.06 


+0.25 
-0.25 


8-9 


n 8fi+0-10 +0.04 +0.06 
'-'■°"-0.10 -0.06 -0.05 


+0.22 
-0.21 


9-10 


n yc+O.ll +0.04 +0.05 
''-'-0.11 -0.06 -0.04 


+0.20 
-0.19 


10-11 


n. Qi +0.16 +0.05 +0.06 
'-'•^-'--0.15 -0.07 -0.05 


+0.25 
-0.25 


11-12 


n Qi +0.19 +0.09 +0.07 
J^-0.17 -0.10 -0.06 


+0.24 
-0.24 


12-13 


n fi8+0-18 +0-05 +0.04 
'-'•'-'°-0.16 -0.07 -0.04 


+0.19 
-0.18 


13-15 


n «Q+0.20 +0.07 +0.04 
'-'•^^-0.18 -0.07 -0.04 


+0.18 
-0.18 



ChiGen Monte Carlo generator [15], which uses a leading-order colour- singlet model, and 
from the NLO NRQCD calculations [3] are shown as hatched bands in Fig. [3] 

Figure |3] also shows the maximum effect of the unknown Xc polarisations on the re- 
sult, shown as the lines surrounding the data points. In the first p"^^ bin, the upper limit 
corresponds to the spin state combination {fnxcii^Xc2) — (0) 2) and the lower limit corre- 
sponds to the spin state combination (1, 1). In all subsequent p'^P' bins, the upper limit 
corresponds to spin state combination (1,2) and the lower limit corresponds to (0,0). 
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p:J^'^ [GeV/c] 



Figure 3: Ratio a{xc2) /^(Xci) in bins of 2<pr^ < 15 GeV/c. The LHCb results, assuming 
the production of unpolarised Xc mesons, are shown with solid black circles and the internal 
error bars correspond to the statistical error; the external error bars include the contribution 
from the systematic uncertainties (apart from the polarisation). The lines surrounding the data 
points show the maximum effect of the unknown Xc polarisations on the result. The upper and 
lower limits correspond to the spin states as described in the text. The CDF data points, at 
y^= 1.96 TeV, are shown with open purple circles [6j. The two hatched bands correspond to 
the ChiGen Monte Carlo generator prediction [15j and NLO NRQCD [3]. 

In summary, the ratio of the (y{Xc2) / <^{Xci) prompt production cross-sections has 
been measured using 36 pb~^ of data collected by LHCb during 2010 at a centre-of-mass 
energy ^/s = 7 TeV. The results are in agreement with the NLO NRQCD model for 
p^^ > 8 GeV/c, whilst there is a discrepancy at low p^^ . Although the ChiGen generator 
describes the shape reasonably well, the data lie consistently above the theory predic- 
tion. This could be explained by important higher order perturbative corrections and/or 
sizeable colour octet terms not included in the calculation. 
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